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ABSTRACT 

Phase-resolved medium resolution VLT spectroscopy of the low mass X-ray binary 
GX 9+9 has revealed narrow C III emission lines that move in phase relative to our 
new estimate of the ephemeris, and show a velocity amplitude of 230±35 km s^ 1 . We 
identify the origin of these lines as coming from the surface of the donor star, thereby 
providing the first estimate of the mass function of /(Mi)^O.22M . Rotational broad- 
ening estimates together with assumptions for the mass donor give 0.07<g^0.35 and 
182^^2^406 km s _1 . Despite a low mass ratio, there is no evidence for a superhump 
in our dataset. Doppler maps of GX 9+9 show the presence of a stream overflow, ei- 
ther in the form of material flowing downward along the accretion disk rim or in a 
similar fashion as occurs in high mass transfer rate cataclysmic variables known as 
the SW Sex stars. Finally we note that the Bowen region in GX 9+9 is dominated by 
CIII instead of NIII emission as has been the case for most other X-ray binaries. 

Key words: accretion, accretion disks - starsundividual (GX 9+9) - X-rays: binaries. 



1 INTRODUCTION 

One of the main aims of optical observations of low mass 
X-ray binaries (LMXBs) has always been to find a signa- 
ture of the donor star, and thereby constrain the masses 
of both components. However, until recently it were mainly 
the LMXBs that show transient outbursts for which it was 
possible to determine the mass function (see e.g. Charles & 
Coe 2006). During a transient outburst the optical emission 
is completely dominated by the X-ray irradiated accretion 
disk, and it is only when the source has returned to quies- 
cence (or at least the X-ray flux has decreased by several 
orders) that the donor star becomes directly visible. That is 
why, until recently, there were no strong constraints on the 
mass of the donor and compact object for most LMXBs that 
are persistently bright, although optical counterparts have 
been known for many systems for over 20 years. 

This situation changed when Steeghs & Casares (2002) 
developed a new technique to detect the companion star in 
persistent LMXBs. Using phase-resolved blue spectroscopy 
they revealed narrow emission line components that were 
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interpreted as coming from the irradiated companion star of 
ScoX-1. This enabled them for the first time to determine 
the mass function. These narrow features were most obvi- 
ous in the Bowen blend (a blend of NIII 4634/4640 A and 
CIII 4647/4650 A lines), which is the result of UV fluores- 
cence due to the hot inner disk for the N III lines or photo- 
ionization and subsequent recombination for the C III lines 
(McClintock et al. 1975). This led to the discovery of such 
signatures in other persistent LMXBs such as X 1822— 371, 
4U 1636-53 and 4U 1735-44 (Casares et al. 2003, 2006), 
or transient sources where the system parameters could not 
be constrained in quiescence such as GX339— 4 or AqlX-1 
(Hynes et al. 2003, Cornelisse et al. 2007). 

GX 9+9 is a persistent LMXB that is bright in both op- 
tical as well as X-rays and is another good candidate for a 
Bowen study. Its compact object is thought to be a neutron 
star due to its X-ray spectral and timing properties (e.g. 
Hasinger & van der Klis 1989), although direct evidence 
such as Type I X-ray bursts or pulsations have not been 
observed thus far. The optical counterpart of GX9+9 was 
identified by Davidsen et al. (1976) as a V—16.6 blue object. 
A 4.2 hour periodic modulation in X-rays was detected that 
was interpreted as the orbital modulation (Hertz & Wood 
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Figure 1. Scarglc normalised power spectrum of the RXTE/ASM 
data of GX 9+9, showing a strong peak at the expected (orbital) 
period of ~4.2 hrs. 



Figure 2. Folded RXTE/ASM X-ray lightcurvc of GX9+9 on 
the orbital period, where phase zero is defined as the minimum 
in the lightcurve. For clarity, the lightcurve is shown twice. 



1988), and confirmed in the optical by Schaefer (1990). From 
its aperiodic variability at very low frequencies, Reig et al. 
(2003) speculated that GX9+9 has an early type (earlier 
than G5) companion. Haswell & Abbott (1994) suggested 
that GX 9+9 is a persistent superhumper, based on its sup- 
posed mass ratio (q<0.28), the optical modulations, and the 
change in the shape of the optical lightcurves. Superhumps 
are periodic optical modulations first observed in SU UMa 
dwarf novae (see Warner 1995), but also in black hole soft 
X-ray transients (O'Donoghue & Charles 1996), and are due 
to an eccentric accretion disk that is undergoing precession 
(e.g. Whitehurst 1988). 

Recent simultaneous X-ray and optical observations by 
Kong et al.(2006) did not detect a significant X-ray/optical 
correlation. Although they could not confirm the X-ray mod- 
ulation at the orbital period, it is present in the RXTE/ASM 
data (Levine et al. 2006a,b). Time resolved X-ray spec- 
troscopy by Kong et al. (2006) showed that a two-component 
spectral model was needed to fit the X-ray emission. The 
best model suggested that an optically thick boundary 
layer which is partly obscured (by either the inner disk or 
thickened disk) and an extended hot corona are present in 
GX9+9. 

In this paper we present VLT medium resolution blue 
spectroscopy of GX9+9, and we derive a new ephemeris 
based on X-ray data. We determine radial velocity curves 
and Doppler maps of the most important lines and use these 
to derive the first constraints on the system parameters. 



2 OBSERVATIONS AND DATA REDUCTION 

On June 24 2003 we obtained two spectra of GX 9+9 in order 
to establish the strength and velocity structure of its Bowen 
blend. This was followed by a phase-resolved study from 25 
to 27 May 2004 using the FORS 2 spectrograph attached to 
the VLT Unit 4 (Yepun Telescope) at Paranal Observatory 
(ESO). During each night in 2004 we observed GX9+9 for 
approximately one full orbit with the 1400V volume-phased 



holographic grism, resulting in a total of 58 spectra (+2 in 
2003) with an integration time of 600 s each. We used a slit 
width of 0.7", giving a wavelength coverage of AA4514-5815 
with a resolution of 70 km s~ (FWHM). We also observed 
GX9+9 once with a wide, 2.0 arcsec slit width. The see- 
ing during these three nights varied between 0.5 and 2.1 
arcsec, and was 1.1 arcsec during the observation with the 
wide slit. The slit was orientated at a position angle of - 
58.5° to include a comparison star in order to correct for 
slit losses. Since we do not have a flux standard, we were 
not able to correct for the instrumental sensitivity curve. 
During the daytime He, Ne, Hg and Cd arc lamp exposures 
were taken for wavelength calibration. We determined the 
pixel-to- wavelength scale using a 4th order polynomial fit to 
20 reference lines giving a dispersion of 0.64 A pixel -1 and 
rms scatter <0.05 A. We de-biased and flat-fielded all the im- 
ages and used optimal extraction techniques to maximise the 
signal-to-noise ratio of the extracted spectra (Home 1986), 
using the routines from the PAMELA package. We also cor- 
rected for any velocity drifts due to instrumental flexure (al- 
ways <5 km s _1 ) by cross-correlating with the sky spectra. 
The subsequent analysis was performed using the package 
MOLLY. 



3 DATA ANALYSIS 
3.1 Orbital period 

The first, and only, optical orbital ephemeris of GX9+9 
was published by Schaefer (1990) over 16 years ago. Unfor- 
tunately, due to the propagation of the phase uncertainty, 
this ephemeris cannot be projected to the present. However, 
Levine et al. (2006a,b) showed that it is possible to detect 
the orbital period in the publicly available X-ray lightcurve 
provided by the All Sky Monitor onboard the RXTE satel- 
lite (RXTE/ASM). Therefore, we decided to derive a new 
ephemeris from this dataset by also searching for the 4.2 
hrs period and determine a new zero point, using the indi- 
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Figure 3. Average normalised spectrum of GX9+9 from our 3 nights of VLT data. We have indicated the most prominent emission 
lines and interstellar diffuse bands (DIB). Furthermore we have indicated a sky line (N/S) that has not been properly removed. 



vidual dwell, full 2-12 keV band-pass observations. We ap- 
plied a barycentre correction to the data using the IDL rou- 
tine HELIO_JD and searched for the orbital period with the 
Lomb-Scargle technique (Scargle 1982). Only one significant 
peak is observed, at 4.19345(2) hrs (see Fig.[T}. This value 
is similar to that of Levine et al. (2006b), but note that 
this period is slightly, but according to the formal errors 
significantly (4.7a), shorter than the photometric period of 
4.1958±0.0005 hrs detected by Kong et al. (2006). In order 
to see if this difference is due to a change in period, we have 
also split the RXTE/ASM dataset in two equal lengths and 
determined the orbital period for each. However, both peri- 
ods are equal, and we conclude that the orbital period has 
not significantly changed over the ~11 yr baseline of the 
RXTE/ASM data. 

We phase-folded all the data on the observed period 
and fitted a sine curve, using a date close to our observations 
as the zero point. We then determined the minimum in our 
folded RXTE/ASM lightcurve to be JD 2453151.464(3), 
and define this as orbital phase zero. In Fig. [2] we show 
the folded lightcurve. Although the lightcurve exhibits a 
modulation of ~4%, comparable to that observed by Hertz 
& Wood (1988), Levine et al. (2006a) reported that the 
amplitude changes over time and even increased up to 18% 
after Jan 2005. Note that we define our orbital phase to be 
half an orbit removed from that used by Schaefer (1990), 
i.e. more in line with the standard definition where phase 
zero corresponds to inferior conjunction of the secondary 
star. This gives a new ephemeris (in HJD) of: 

T=2453151.464(3)+0.1747272(8)E 

which we will use in the remainder of this paper. In Sects. 3.3 
and 3.5 we will determine the ephemeris in three different, 
independent, ways. They are all consistent with the above 
ephemeris, giving us confidence that it represents a reliable 
orbital phasing. However, we must always be careful with the 
assumption that the minimum in the lightcurve corresponds 
to inferior conjunction of the secondary star. 



3.2 Flux calibration and spectrum 

In order to correct for slit losses due to the variable seeing 
we used a comparison star placed on the slit. We divided 
all comparison star spectra by its wide-slit spectrum, and 
fitted them with a low-order spline. The spectra of GX 9+9 
were subsequently divided by the comparison star splines 
(corresponding to the same observation) in order to get the 
final spectra corrected for slit losses. Note that since we do 
not have a flux standard all our spectra have relative fluxes. 

Figure|3] shows the average spectrum of GX 9+9 nor- 
malised to unity. It is dominated by high excitation emis- 
sion lines from Hell A4686, A5411 and Bowen AA4630-4650. 
A strong absorption feature with emission structure super- 
posed is present around H/3. Several absorption lines can 
also be observed that are mainly due to interstellar bands 
(indicated in Fig.[3]as DIB). A close inspection of the Bowen 
region shows that there is little emission below 4640 A, i.e. 
the part of the blend where all the dominant N III fluores- 
cence lines should be. In order to quantify this absence of 
NIB we measured the centroid of the Bowen region using 
a Gaussian, and find a central wavelength of 4648.0+0.2 
A (while Hell has a central wavelength of 4686.7+ 0.2 A). 
This central wavelength is much higher than was the case for 
e.g. V801 Ara (4643.2 A) and V926Sco (4641.5 A; Casares 
et al. 2006). Although there might be some NIII present, 
it does suggest that the Bowen blend is dominated by CHI 
in GX9+9. We therefore checked where other strong CIII 
transitions should occur in our spectra. There should be one 
emission line at 4860 A, with a transition strength compa- 
rable to those in the Bowen region. Although this is close 
to H/3 we will show in Sect. 3.5 that we have most likely de- 
tected it, although the presence of Hell A4859 can also not 
be ruled out. Another strong CIII transition occurs at 5305 
A, and is also present in our spectrum (indicated in Fig. [3}, 
although the width of the line suggests that it is a blend 
with other unknown lines. Finally, we checked if any strong 
NIII or He I lines are present in our spectral range (apart 
from the NIII lines in the Bowen region). Although there 
should be 3 strong N III lines at 4867, 5059 and 5320 A and 
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Figure 4. Average lightcurves of GX9+9 with the maximum 
normalised to unity (also for H/3, despite negative flux levels) and 
folded on the orbital period of 4.19 hrs for the continuum, H/3, 
Hell A4686 and the Bowen blend respectively. For clarity each 
lightcurvc is plotted twice. 



2 He I lines at 4922 and 5015 A, none of these appear to be 
present. 



3.3 Optical lightcurves 

We used the orbital ephemeris derived above (Sect. 3.1) to 
compute the orbital phase of each spectrum. We integrated 
the continuum flux from each individual spectrum, phase- 
folded them into 20 bins, and normalised the peak flux to 
unity in order to create the lightcurve that is shown at 
the top of Fig.[4] The first thing to note is that the min- 
imum in the optical continuum lightcurve corresponds very 
well with orbital phase as defined by the X-ray modula- 
tion, and thereby strengthening our interpretation that the 
X-ray ephemeris represents inferior conjunction of the sec- 
ondary. Furthermore, we also note that the flux varies by 
~20%, comparable to previous observations (e.g. Schaefer 
1990, Kong et al. 2006). 

We have also subtracted the continuum lightcurve from 
each spectrum in order to derive light curves from the main 
emission or absorption lines, and these are also shown in 
Fig. [4] We do note that on average the flux in the lines vary 
by 40%, i.e. the variation is twice as strong as in the contin- 
uum. The flux of the H/3 absorption line shows more or less a 
sinusoidal variation, with minimum (i.e. maximum absorp- 
tion) occurring around orbital phase 0.5 (in anti-phase with 
the continuum emission) . The variation of He II A4686 and 
the Bowen blend is more complex. Hell A4686 appears to 
be a combination of the H/3 and continuum lightcurves, i.e. 
showing a minimum around orbital phase 0.5 and a second 
minimum at phase 1.0. The Bowen emission, on the other 
hand, looks similar to the continuum emission, although its 
single minimum is at a slightly later orbital phase of ~ 1.1. 
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Figure 5. Trail of the Bowen blend and Hell A4686 line (top 
panel), and H/3 (bottom panel), folded on the orbital period and 
plotted twice for clarity. The numbers indicated in the trails cor- 
respond to the component numbers in Table[T] The light areas in 
the H/3 trail are the absorption feature, while dark regions are the 
emission features. 



3.4 Trailed spectra and radial velocity curves 

We have phase-folded all data into 20 bins, and plotted a 
trail of the region around He II A4686 and the Bowen blend 
in Fig. [5] (top panel). The HeIIA4686 line clearly consists of 
two separate components that we have indicated with 1 and 
2 in Fig. [S] and correspond to the component number in Ta- 
ble[T] Both components show a sine-like modulation in anti- 
phase with each other. Note that these oppositely phased S- 
waves also produce double peaks, but these are completely 
different from the more conventional disk dominated line 
profiles where the peaks move together in phase. There is 
no strong broad component visible, indicating that most of 
the He II A4686 observed is coming from two distinct regions 
in the binary. We estimated the centre of each emission com- 
ponent as a function of orbital phase using two Gaussians, 
and then fitted each component with a sine curve, obtaining 
the parameters given in Table[T] Fig. [5] already suggests that 
component 2 has twice the velocity amplitude of component 
1, and this is confirmed in Table[T] 

We also examined the region around H/3 in more detail, 
and show a trail in Fig. [5] (bottom panel). Due to the pres- 
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Table 1. Kinematics of GX9+9 emission line structures obtained 
by best fitting sine-curves to the double peaked structure observed 
in He If A4686 and H/3 regions. The component numbers for He II 
correspond to the numbers indicated in the Hell trail of Fig. [5] 
and H/3 component 2 has the same origin as component 2 in He II. 





component 


K 


7 






# 


(km s _1 ) 


(km s — 1 ) 




Hell 


1 


217±5 


97±5 


0.07±0.05 


Hell 


2 


402±3 


91±3 


0.55±0.05 


cm 


3 


260±5 


87±5 


0.02±0.05 


H/3 


2 


398±8 


100±6 


0.49±0.05 



K =semi-amplitude,7= off-set to rest-wavelength, 
c/>=phase zero 



ence of the absorption feature this region is more complex. 
However, we do notice that the emission feature also appears 
to consist of two components (indicated with 2 and 3), more 
or less in anti-phase with each other. These components are 
clearly visible in the individual spectra (see e.g. Fig. [6}, and 
made us confident that we could again fit two Gaussians to 
the emission components (see also Table[T]for the best fit- 
ting sine curves, which we have labelled CIII and H/3). We 
note that the strongest feature (called component 2) appears 
to originate from the same region as component 2 of Hell 
A4686, and we will discuss this in more detail in Sect. 3.5. 
On the other hand, if we assume that component 3 is also 
H/3 emission, its average velocity, j, is completely different 
compared to all the other components listed in Table[T] (7 
would be -16 km s _1 compared to ~95 km s _1 ). However, if 
we assume this line is CIII A4859.65 and/or Hell A4859.3 
(see Sect. 3.2) 7 becomes comparable with that derived for 
the other emission lines. We, therefore, conclude that com- 
ponent 3 is most likely tracing a combination of C III and 
He II. In Sect. 3.5 we will show that this region is most likely 
dominated by C III and H/3 emission, and despite its com- 
plicated appearance this trail is very similar to Hell A4686. 

In the H/3 trail (Fig. [5} it appears that no net emission 
components are present at orbital phase 0.5. However, given 
that the absorption feature hardly moves with orbital phase, 
we expect the emission feature to be at the deepest point of 
the absorption feature at orbital phase 0.5. If the peak flux 
of the emission feature is below the continuum level it would 
be difficult to detect its signature in the trail. In Fig. [6] we 
show the average spectrum around orbital phase 0.5, and it 
is clear that the H/3 and CIII/Hell emission components 
are still present (but have peak fluxes below the continuum 
level) . 

We also created trails of each individual night by folding 
them on the ephemeris determined in Sect. 3.1. If GX9+9 
is a superhumper, it would most likely show a similar spec- 
troscopic signature as other superhumping LMXBs, namely 
a movement in the nightly mean emission line profiles (see 
e.g. Zurita et al. 2002). Therefore, we did test if movement 
is present, but no shift in the line centroid was detected be- 
tween the different nights. Also a comparison between the 
2003 and 2004 spectra shows no clear shift that would un- 
ambiguously point toward a spectroscopic signature of a su- 
perhump in GX9+9. 

In order to determine the radial velocity of the compact 
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Figure 6. Average spectrum of the H/3 region of GX9+9 around 
orbital phase 0.5, clearly showing that the main emission compo- 
nents superposed on the broad absorption are still present. 



object one can apply the double- Gaussian technique (Schnei- 
der & Young 1980). We must note that this technique only 
works for a standard double peaked emission profile disk, 
which is clearly not observed in GX9+9 (see Fig.[5jl. We 
still decided to apply the technique to He II A4686 using ve- 
locities that are as far as possible from the prominent emis- 
sion features in the trail, in the hope that their contribu- 
tion is negligible. We used a double Gaussian bandpass with 
FWHM of 100 km s _1 and Gaussian separations between 
600 and 1200 km s - . Fig. [7] shows that as we move away 
from the line core, the K\ velocity exponentially declines to 
120±10 km s~ , while the systemic velocity, 7, moves be- 
tween 100 and 123 km s _1 . The phasing for the compact 
object showed a small change over Gaussian separation, but 
corresponds well with the value predicted by the ephemeris 
derived above (i.e. ~0.5). Although we expect this kind of 
behaviour for the K\ velocity we are still hesitant to iden- 
tify this as the true motion of the compact object, and we 
should treat the result with caution. 



3.5 Doppler maps 

Doppler tomography is a powerful tool for probing the struc- 
ture of the accretion disk (Marsh & Home 1988). It makes 
use of all the available data and is effective for emission fea- 
tures that are too weak to be identified in the individual 
spectra, by resolving the distribution of the line emission in 
the co-rotating frame of the binary system. The definition 
of orbital phase in Doppler maps corresponds to inferior 
conjunction of the secondary (Marsh 2001), and is one of the 
main reasons why we changed the definition of our phase 
compared to that adopted by Schaefer (1990). This should 
project a signature of the donor star along the positive y-axis 
(Steeghs & Casares 2002). 

Fig. [8] shows the Doppler maps of the three most im- 
portant emission line regions, i.e. Hell A4686, the Bowen 
blend and H/3. Since there are no other emission lines present 
around Hell, we could use the standard reduction tech- 
niques with a systemic velocity of 90 km s . As expected 
from the trails (see Sect. 3.4) two main emission regions can 
be observed (Fig. [8j/top- left), more or less in anti-phase with 
each other. The bright region in the top half corresponds 
with the position where the donor star and accretion stream 
are expected to be located. However, the brightest feature 
is the extended region in the bottom-left quadrant of the 



6 R. Cornelisse et al. 



a (kn 



Figure 7. Diagnostic diagram for GX 9+9 of He II A4686 to follow 
the estimates of the radial velocity of the compact object, K\, and 
the systemic velocity, 7, as a function of Gaussian separation, a. 



Doppler map, a position where no enhanced emission is ex- 
pected for a 'normal' X-ray binary. 

Using the correct systemic velocity should produce the 
sharpest defined feature in the top part of the Doppler map. 
We therefore tried different velocities with steps of 10 km s _1 
around 90 km s _1 . However, in line with Sect. 3.4, we found 
that 90 km s _1 gives the sharpest result for the compact 
spot in the top part of the map, and again interpreted this 
as the systemic velocity. Furthermore, we created Doppler 
maps of Hell A4686 for each individual night in order to 
explore any significant change. All maps looked identical, 
and no long-term change in the emission line profile appears 
to be present. 

The Doppler map of the Bowen region (Fig. [5J/top- right) 
was created using standard reduction techniques, but fitting 
the three strongest C III lines (4647,4650 and 4665 A) simul- 
taneously using relative strengths as indicated by McClin- 
tock et al. (1975). The most pronounced feature in the map 
is a sharp spot in the top half of the map, at a position 
similar to that in the He II map and where the donor star is 
expected. Fitting a 2-dimensional Gaussian to this spot we 
measure a velocity of 250±25 km s _1 . This 25 km s" 1 is the 
formal 1-sigma fitting error in the centroid position of our 2- 
dimensional Gaussian. Although the spot is offset from the 
x-axis by 35 km s _1 , this only corresponds to a phase offset 
of 0.02 which is smaller than the error in the phase origin 
(0.04). The spot location is thus consistent with being along 
the axis connecting the neutron star with the donor star. 
We also rotated the map by 0.02 orbital phase in order to 
align the spot with the x-axis. Since we measured a similar 
velocity as above, we concluded that the error introduced 
by the off-set is negligible. 

The most complex region is H/3, where C III A4860/He II 
A4859 and an H/3 emission line are superposed on a deep 
absorption feature. Since it is currently not possible to di- 
rectly create Doppler maps with negative pixel values we 
used a non-standard reduction technique. First we created 
a Doppler image that consisted of a Gaussian with a peak 
and width comparable to the observed absorption feature. 
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Figure 9. Average spectrum of the Bowen (top) and H/3 (bottom) 
regions of GX9+9 in the rest-frame of the donor star. We have 
indicated the dominant contributing lines in these regions. 



From this image we created artificial data that has the same 
properties as the observations, i.e. creating a trail of Gaus- 
sians. We then added the artificial data to the observations 
in order to remove the absorption feature. We then used 
standard analysis techniques to create Doppler images for 
both the C Ill/He II and H/3 simultaneously. Finally, we sub- 
tracted the artificial map of the Gaussian to produce the 
maps shown in Fig. [8] (bottom) . Although we were able to 
more or less disentangle the different components, we must 
be cautious about the details of the two maps. 

The CIII/He II Doppler map (Fig.EJ/bottom-right) 
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Figure 8. Doppler maps of Hell A4686 (top-left), the Bowen region (top-right). The bottom panel shows the Doppler deconvolution of 
the H/3 (bottom-left) and C III A4860 components as explained in Sect. 3.5. In Hell, Bowen and CIII the white parts indicate the zero 
level of the flux, while in H/3 the white parts indicate regions of most absorption. Indicated on the Bowen and C III maps are the Roche 
lobe, gas stream leaving the LI point, and the Keplerian velocity along the stream for the set of parameters q=0.l6 and K2=264 km 
s — 1 discussed in Sect. 4.2. 



looks remarkably similar to that of the Bowen region, sug- 
gesting that the line is indeed dominated by CIII emission. 
The main feature is a single localised spot in the top half 
of the map, that is at a velocity of 210±25 km s _1 (and an 
off-set from the x-axis of ~55 km s" 1 , still within our phase 
zero uncertainty), consistent with the spot in the Bowen 
map. The H/3 map (Fig.[8j/bottom-left) looks more com- 
plex. All absorption is attributed to H/3 and is present at 
all orbital phases, although it is most prominent in the left 
part of the map around y—0, £=-500 km s — \ Note that 
in this map the lightest regions indicate negative flux lev- 
els, in sharp contrast with the other maps where lightest 
regions represent the continuum level. However, the most 
prominent feature in the H/3 map is created by the emission 
lines, namely the extended region in the bottom half of the 
map that has a comparable shape and position to that in 



the He II map. However, this similarity between the H/3 and 
Hell maps is rather curious, since in other LMXBs these 
lines can show rather different characteristics (e.g. Steeghs 
& Casares 2002). Although this might be explained by con- 
tamination of He II A4859, we would have expected that this 
would mainly affect the C III map. We must, therefore, con- 
clude that this emission region in the bottom half of the H/3 
map is likely real. 



4 DISCUSSION 

We have presented medium resolution optical spectroscopy 
of GX9+9. One of the most remarkable characteristics of 
GX9+9 is the strong CIII emission in its spectrum, while 
NTH is hardly detected. To our knowledge, no other X-ray 
binary is known that shows such an unusual CIII/NIII ra- 
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Table 2. Derived system parameters for GX9+9. 



Parameter 




Parameter 




P orb (hrs) 


4.1934528(20) 


f(M) (M a ) 


(^)0.22±0.03 


T (HJD) 


2 453 151.464(3) 


q (M 2 /Mi) 


0.07-0.35 


7 (km s — 1 ) 


90±7 


K 2 (km s" 1 ) 


182-406 


Kern (km s" 1 ) 


230±35 







tio in the Bowen region. For example, in X1822— 371 and 
ScoX-1 the narrow Carbon and Nitrogen lines in the Bowen 
region were more or less of comparable strength, and no 
other strong CIII/NIII lines could be detected (Casares 
et al. 2003, Steeghs & Casares 2002). This could give us im- 
portant constraints on the evolutionary status of GX9+9. 
However, we can exclude a more evolved core in the donor 
star of GX9+9, since this would lead to an overabundance 
of NTH due to processing by the CNO cycle as has been 
observed in e.g. XTE Jl 118+480 Haswell et al. (2002). The 
only other systems that show prominent Carbon emission 
lines in their spectra are some of the ultra-compact X-ray 
binaries (Nelemans et al. 2006). However, in these systems 
the donor star is thought to be a white dwarf that is defi- 
cient of both H and He. These systems are clearly different 
to GX 9+9, which is not an ultra-compact X-ray binary and 
has both H and He in its spectrum. 



4.1 Orbital Period and Long Term Variations 

Following Levine et al. (2006a, b) we have updated the 
ephemeris of GX9+9 using the RXTE/ASM data. Given 
the coincidence of optical minimum with the minimum in the 
X-ray lightcurve, as well as the presence of Bowen emission 
spots on the y-axis consistent with the donor star, we feel 
confident that this ephemeris tracks the correct orbital phase 
of the binary. Levine et al. (2006a) showed that the strength 
of the modulation in GX 9+9 increased from 6% before Jan 
2005 to 18% after this date (and until at least June 2006). A 
possible explanation for this change in modulation is given 
by the folded RXTE/ASM lightcurve shown in Fig.0 There 
appears to be a partial eclipse around orbital phase 0, i.e. 
at inferior conjunction of the secondary. The most likely re- 
gion that could be eclipsed by the secondary is the extended 
accretion disk corona that Kong et al. (2006) observed in 
GX 9+9, suggesting that the corona slowly changes size. One 
option could be that the orbital modulation in X-rays dis- 
appears when the corona becomes too small to be eclipsed 
by the secondary. An alternative option is that the corona 
is always eclipsed, but when it increases in size its fractional 
contribution to the total X-ray flux increases thereby lower- 
ing the depth of the eclipse. If the size of this corona could 
be determined, we would be able to distinguish between the 
different scenarios, and could also constrain the inclination 
of GX9+9. However, the fact that this partial eclipse occurs 
suggests a reasonably high inclination. 

Haswell et al. (2001) proposed that GX9+9 is a persis- 
tent superhumper. Superhumps are periodic optical modula- 
tions first observed during outburst of SU UMa dwarf novae 
(see Warner 1995 for a review), and their most important 
property is that the photometric period is a few percent (1- 




Figure 10. Constraints on q and K2 for GX9+9. K2 must be 
larger than K C m, and a disc opening angle of a=0 gives the upper- 
limit for K2 ■ The rotational broadening derived from the emission 
line widths sets a lower limit to q, while the largest possible main 
sequence star that fits the Roche-lobe of GX 9+9, combined with 
its maximum allowed broadening sets the upper-limit to q. All 
curves indicate 95% confidence levels. 



7%) longer than the orbital period. This is thought to be due 
to prograde precession of an eccentric accretion disk, giving 
rise to slightly longer periods (with respect to the orbital pe- 
riod) for the accretion disk. The main criterion for creating 
an eccentric disk is an extreme mass ratio (M2/Mi<0.33). 
Comparing our X-ray period with the photometric period 
derived by Kong et al. (2006) we note that the photomet- 
ric one is 0.0564 hrs longer, corresponding to only ~1% of 
the orbital period. Although this difference is small, it is 
still significant according to the formal errors (4.7<r), and 
we therefore further investigated the possibility of a super- 
hump in GX9+9. 

If GX 9+9 is a superhumper, we would also expect sim- 
ilar spectroscopic characteristics for GX9+9 to that of the 
known superhumping LMXB XTE J1118+480, namely the 
movement of the nightly mean emission line profiles (e.g. Zu- 
rita et al. 2002). Such movement is not observed in GX9+9, 
arguing against a superhump interpretation. Furthermore, 
also the so-called e(q)-q relation, where q is the mass ratio 
and e(q) the fractional period excess of the superhump over 
the orbital period, suggest that the photometric period can- 
not be due to a superhump. If the difference in photometric 
and X-ray period is real, the empirical relation by Patterson 
et al. (2005) suggests that g~0.003, an unrealistically small 
value (as we will show in Sect. 4.2). This could indicate that 
the errors on the orbital periods have been underestimated, 
and that the X-ray and photometric period are similar. On 
the other hand we cannot completely exclude the possibility 
that the movement of the nightly mean emission line pro- 
file is hidden by the strong emission feature in the lower-left 
quadrant of the Doppler maps (Fig. [5}, and that the e(q)- 
q does not apply to persistently bright LMXBs. The only 
thing we therefore conclude is that there are no indications 
of a superhump in our dataset, and that a better determi- 
nation of the photometric period is needed to find out if the 
difference with the X-ray period is real. 
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4.2 Donor Star and System Parameters of 
GX9+9 

Thus far narrow lines in the Bowen region have been de- 
tected in ScoX-1, X 1822-371, GX339-4, 4U 1636-53, 
4U 1735-44 and AqlX-1 (Steeghs & Casares 2002; Casares 
et al. 2003, Hynes et al. 2003, Casares et al. 2006, Cor- 
nelisse et al. 2007). We can now also add GX9+9 to this 
list. Phase zero of the radial velocity curve of these narrow 
lines in GX9+9 correspond well with the minima in the 
X-ray and optical lightcurves, suggesting that these lines 
arise close to or are connected with the donor star. It has 
been proposed for the other sources that these narrow lines 
arise on the surface of the secondary, and especially for the 
eclipsing LMXB X 1822—371 this connection could unam- 
biguously be claimed (Casares. 2003). Since it is possible, 
for a valid set of system parameters (see below), to draw a 
Roche lobe on the Bowen and C III Doppler maps that en- 
compasses most of the observed narrow spots (see Fig. [5}, we 
propose that also in the case of GX 9+9 these lines arise on 
the donor star. If our suggestion is correct, the narrow fea- 
ture in the Doppler maps of both the Bowen region and C III 
A4860 (see Fig. [8]) must be produced somewhere on the inner 
hemisphere of the companion, and has a Kem of 230±35 km 
s _1 (the average value derived from the compact spots in 
the Doppler maps of Hell A4686, Bowen and CIII A4860). 
However, given that only the irradiated side facing the com- 
pact object is contributing to this emission, this should be 
considered a lower-limit to the velocity semi-amplitude of 
the centre of mass of the donor Ki- 

In order to derive some constraints on the system 
parameters of GX9+9, we can start by determining a 
lower limit on the mass function of f(M)=M 1 sin 3 i/(l + 
g) 2 >0.22±0.03M Q (where q=M 2 /M 1 is the mass ratio and 
i the inclination of GX9+9). This lower limit corresponds 
with the assumption that K eln =K2, i.e. that the narrow lines 
are produced on the poles of the donor star. In order to 
find the real K2 we need to determine the A"-correction, 
which is dependent on the disk flaring angle (a), q, and 
weakly on i (see Munoz-Darias et al. 2005). Using an ir- 
radiation binary code, Munoz-Darias et al. determined 4th 
order polynomials for different values of a and i to estimate 
the Tf-correction. Unfortunately we do not have any solid 
estimates of Ki, since the double Gaussian method used 
in Sect. 3.4 is not valid for an accretion disk that is domi- 
nated by a hotspot. Therefore we cannot directly apply the 
_R"-correction, and must therefore resort to several assump- 
tion in order to further constrain the system parameters. 
Although we will point out what the major source of uncer- 
tainty is for each of our assumptions, we do think that they 
are currently the best we can do. 

First of all, the simplest assumption is that K2^K em , 
and show the 95% confidence limit in Fig. [TO] Furthermore, 
we can use the polynomials by Munoz-Darias et al. (2005) 
for a=0° (and i=40°) to determine the upper- limit to K2 
as a function of q. Again we show the 95% confidence curve 
in Fig.fTUl 

We can determine a lower-limit to q by estimating 
a lower limit to the rotational broadening, Usini, of the 
emission lines produced on the surface of the donor (Wade 
& Home 1988). However, this does mean that we assume 
that all of the line broadening is due to rotation, which 



o .0 




Velocity (km/s) 

Figure 11. He II A4686 Doppler map showing the Roche lobe of 
the secondary, the gas stream leaving the LI point, and the Kc- 
plerian velocity along the stream for two different sets of system 
parameters. The circle indicates the Keplerian velocity at the edge 
of the disk. Top: System parameters derived for the assumption 
that K\ derived with the double Gaussian technique is correct. 
Bottom: System parameters derived using the assumption that 
the emission feature in the bottom left quadrant is part of the 
accretion disk. 



might not be the case. We created an average spectrum in 
the rest- frame of the secondary. Fig. [9] shows the two most 
prominent regions (around He II A4686 and H/3) of the rest- 
frame spectrum. Unfortunately, most lines are either too 
weak, are located in very complicated regions (e.g. Bowen 
and H/3 regions) or have a component that does not orig- 
inate on the secondary (e.g. Hell A4686). We decided to 
estimate the FWHM of the two CIII lines in the Bowen 
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region, C III A4859, and also the H/3 absorption line. All 
lines give a similar FWHM that is on average 153±26 km 
s _1 . However, we do note that this includes the effect of 
the intrinsic broadening due to the instrumental resolution 
of 70 km s _1 . To take this into account we broadened a 
strong line in our arc spectrum using a Gray rotational pro- 
file (Gray 1992) until we reached the observed FWHM. We 
assumed no limb-darkening since the fluorescence lines oc- 
cur in optically thin conditions. We found that a rotational 
broadening of 102±15 km s _1 reproduces our results. Using 
Kotsim=0.462if 2 g 1/3 (1 + q) 2/3 (Wade & Home 1988), and 
Vsini^73 km s _1 (95% confidence) we can draw this curve 
in Fig.fTOl 

We also show the limit corresponding to the assump- 
tion that the donor star is a zero-age main sequence star 
that fits in a 4.2 hrs period Roche lobe. Although there is 
no guarantee that the donor star in GX9+9 satisfies that 
assumption (see e.g. Schenker & King 2002 for alternative 
scenarios), we include it in Fig. [10] as a likely constraint on 
the system parameters (Casares et al. 2006). This leads to a 
donor star mass ^0.47Mq and together with Mns^1-35Mq 
gives g^0.35. Note that this implies that the lower limit 
on K2 derived from the compact spot in the Bowen and 
C HI Doppler maps is always below the rotational broaden- 
ing curve determined from the widths of the emission lines. 

By estimating the upper-limit to the rotational broad- 
ening for a O.47M0 main sequence star donor we derive the 
final limit shown in Fig.QJJ From Paczynski (1974, 1983) we 
can determine an upper-limit to the inclination of ^65°, i.e. 
where the donor star would just eclipse the compact object. 
Since Kong et al. (2006) did not detect any orbital variation 
in X-rays, and their X-ray spectral modelling showed that 
the emission coming from the inner accretion disk or bound- 
ary layer is a significant fraction of the total flux (~25%), we 
think our assumption on the inclination is valid. This gives 
a maximum rotational broadening of l/sini=123 km s~ , as 
indicated in Fig. 1101 

The values for the system parameters of GX9+9 we 
have derived are the best limits possible, and these are shown 
in Fig. [10] and tabulated in Table[2] In the remainder of this 
section we will speculate in two different ways about further 
constraints on the system parameters using more uncertain 
assumptions. Although they lead to very different system 
parameters, we will use the results of these speculations for 
our discussion of the accretion disk structure in Sect. 4.3. 

Although the double Gaussian technique is not an ap- 
propriate way to determine K\ in the case of GX9+9, the 
value derived is within the acceptable limits. We can there- 
fore speculate that this is the correct K\ and see how this 
limits the system parameters. Again using the polynomi- 
als by Munoz-Darias (2006) we can determine the largest 
JT-correction, which leads to K2 ^369 km s _1 and thereby 
g^0.32. Since we already established that q must be smaller 
than 0.35 (see Fig.QJjJ, we derive a lower limit to K3 of 343 
km s~ x in this case. We note that these system parameters 
are all toward the upper-limits of the allowed area in Fig. 1101 
We use the lower limits of our derived parameters (q=0.35 
and ^2 =343 km s" 1 ) to draw the Roche- lobe of the sec- 
ondary, stream trajectories and outer edge of the accretion 
disk on the Doppler map of Hell and show this in Fig. 1111 
(top panel). We have indicated the Roche lobe of the sec- 
ondary and the outer edge of the accretion disk as solid lines. 



The crosses indicate the velocities of the centre of mass, and 
the two stars, while the dotted lines indicate the ballistic 
trajectory of the material leaving the LI point and the Kep- 
lerian velocity the ballistic stream passes through. It is clear 
from this figure that in this case the emission feature in the 
lower-left quadrant is at too low velocities to be associated 
with the accretion disk. 

Another way to estimate the system parameters is if we 
speculate that the emission feature in the lower-left quadrant 
must be part of the accretion disk. Furthermore, we can also 
assume that the feature close to the donor star signature 
is part of the accretion stream. These constraints imply a 
small K2 (to make the velocities of the outer edge of the 
disk similar to the feature in the lower-left) and also a small 
q (the Roche-lobe of the secondary must be small to also 
let the stream pass through the feature close to the donor 
star). Assuming that our estimate of Vsmi is correct, there is 
a possible combination of q and K2 at the lower limits of our 
allowed region (g=0.16 and 7^2=264 km s" 1 for a=12° and 
z=40°) that fulfils all our criteria, and we show the result in 
Fig.HU (bottom). 

Note that we do not claim that either of the two sets of 
system parameters shown in Fig. [TJJ must be correct. Given 
the large area of possible solutions in Fig.[TU] it would be 
very unlikely that they are, but they nicely illustrate that 
using different assumptions, both not too unreasonable, it is 
possible to derive completely different system parameters. In 
order to really constrain the system parameters of GX9+9, 
it is necessary to unambiguously determine K\, as was for 
example possible for X 1822—371 or V801 Ara (Jonker & van 
der Klis 2001; Casares et al. 2006). Unfortunately, this is go- 
ing to be difficult for GX9+9, since thus far it has not been 
possible to directly observe a signal coming from the com- 
pact object (in the form of pulsations or Type I X-ray bursts) 
or even the most inner regions of the accretion disk (in the 
form of high frequency quasi-periodic oscillations (e.g. van 
der Klis 2006). 

4.3 Accretion Disk Structure 

GX 9+9 shows several interesting characteristics that are not 
commonly observed in LMXBs. Most important of all is the 
presence of a strong emission feature in the lower-left quad- 
rant of both the Hell A4686 and H/3 Doppler maps (see 
Fig. [5}, while another uncommon feature is the extremely 
broad absorption feature around H/3 that is present at all 
phases. 

In many LMXBs, the Hell Doppler map displays a 
classic disk-like structure, with sometimes a spot in the 
upper-left quadrant due to the accretion stream impacting 
the disk. However, having emission predominantly in the 
lower-left quadrant of the Doppler map is not unique for 
GX9+9. There are at least two LMXBs that show a similar 
Hell Doppler map, the transient sources XTE J2123— 058S 
(Hynes et al. 2001) and GRO J0422+32 (Casares et al. 
1995). In particular, a comparison between the Hell A4686 
trails of GX9+9 and XTE J2123-058 shows that they are 
strikingly similar to each other. Hynes et al. (2001) com- 
pared the properties of XTEJ2123-058 with that of the 
SW Sex objects, a subclass of nova-like cataclysmic variables 
(see e.g. Thorstensen et al 1991). Although there are several 
competing models, it is generally thought that SWSex ob- 
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jects have some sort of stream (that is not connected to the 
accretion disk), that is either overflowing the disk or being 
propelled away, that gives rise to their peculiar characteris- 
tics (see e.g. Hoard et al. 2003, Hellier 2000, Home 1999). 

In Fig. [11] we show various system characteristics su- 
perposed on the Hell Doppler map for two different sets of 
system parameters. The important thing to notice in Fig. 1111 
(bottom) is that there are combinations of valid system pa- 
rameters possible where the emission feature in the lower- 
left quadrant overlaps with the expected velocities for the 
outer edge of the accretion (although some emission it is still 
at sub-Keplerian velocities). This suggests that most of the 
He II emission is produced by material that is carried down- 
stream after the stream impacts the accretion disk. Such 
behaviour was also observed in EXO 0748—676 where es- 
pecially emission from the higher ionization emission lines 
showed up further downstream (Pearson et al. 2006), at a 
position comparable to the Hell emission in GX9+9. How- 
ever, there are also many combinations of the system pa- 
rameters possible where all Hell emission in GX9+9 is at 
sub-Keplerian velocities (as illustrated in Fig. llip . and that 
disk overflow must take place. In this case it would be very 
similar to XTE J2123— 058, and it would then be the second 
LMXB to show SWSex behaviour (Hynes et al. 2001). Un- 
fortunately, currently we cannot distinguish between both 
scenarios, and must conclude that both are just as likely. 

The presence of the broad absorption component 
around H/3 in GX9+9 is another curious feature. Although 
other LMXBs also show such a feature, e.g. X 1822—371, 
SerX-1 and XTE J2123— 058 amongst others (Casares et al. 
2003, Hynes et al. 2004, Hynes et al. 2001), it is usually 
not as strong as observed in GX9+9. An exception is the 
X-ray transient GRO J0422+32 (that showed a similar He II 
A4686 Doppler map as GX9+9) where the absorption fea- 
ture around H/3 (and H7) shows a similar strength (Casares 
et al. 1995). Although this might be a coincidence, since 
XTEJ2123— 058 does not show such prominent absorption 
(Hynes et al. 2001), it is curious. An explanation could be 
a high inclination for GX9+9 that causes absorption by 
the optically thick parts of the inner accretion disk. This 
is supported by the shape of the folded X-ray lightcurve 
(Fig. [21 and the fact that Kong et al. (2006) concluded that 
a significant fraction of the inner accretion disk or boundary 
layer is not visible due to obscuration by the inner or thick- 
ened disk. Interestingly, XTE J2123-058 and many SWSex 
objects also have a high inclination (Casares et al. 1998, 
Knigge et al. 2000). Although a high inclination has also 
been suggested for GRO J0422+32 (Kato et al. 1995), there 
are indications that it is lower (e.g. Gelino & Harrison. 2003; 
Reynolds et al. 2007). 



5 CONCLUSIONS 

We have presented a spectroscopic dataset of GX9+9, and 
using the Bowen fluorescence technique have for the first 
time detected narrow lines that we interpret as coming from 
the donor star. If true, this allowed us to determine a lower 
limit on the mass function of the system of 0.22Mq and 
182^A"2^406 km s . In order to constrain the system pa- 
rameters more information is necessary, in particular on K\ . 
However, given that thus far no unambiguous detection of 



the compact object has been made this is going to be diffi- 
cult. 

We have updated the ephemeris of GX9+9 using 
RXTE/ASM X-ray data, and the folded lightcurve suggests 
a partial eclipse of the extended accretion disk corona by the 
secondary. The Doppler maps of GX 9+9 show the presence 
of stream overflow. This is either of material that is carried 
downward along the rim of the accretion disk, or overflow- 
ing material that is at sub-Keplerian velocities comparable 
to XTE J2123-058 and SWSex systems (Hynes et al. 2001). 
Finally, GX9+9 shows a curious CIII/NIII ratio compared 
to other LMXBs, which could provide an important hint 
about the evolutionary status of this system. 
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